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ABSTRACT 

We present a kinematic study of the evolution of coronal mass ejections 
(CMEs) in the solar wind. Specifically, we consider the effects of: (1) spheri- 
cal expansion; and (2) uniform expansion due to pressure gradients between the 
Interplanetary CME (ICME) and the ambient solar wind. We compare these 
results with an MHD model, which allows us to isolate these effects from the 
combined kinematic and dynamical effects, which are included in MHD models. 

They also provide compelling evidence that the fundamental cross section of so- 
called “force-free” flux ropes (or magnetic clouds) is neither circular or elliptical, 
but rather a convex-outward, “pancake” shape. We apply a force-free fitting to 
the magnetic vectors from the MHD simulation to assess how the distortion of 
the flux rope affects the fitting. In spite of these limitations, force-free fittings, 
which are straightforward to apply, do provide an important description of a 
number of parameters, including the radial dimension, orientation and chirality 
of the ICME. 

Subject headings: Sun: coronal mass ejections (CMEs)- Sun: activity-Sun: 
corona-Sun: magnetic fields-solar wind 


1. Introduction 

Magnetic clouds (MCs) or flux ropes (we will not distinguish between the two terms 
here) are transient solar wind structures characterized by: (1) strong magnetic fields; (2) large 
coherent rotations in the magnetic field vector; and (3) low proton temperature (Burlaga 
et al. 1981). Since their discovery in solar wind data more than 20 years ago, they have 
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been a major focus of heliospheric research. The commonly held view is that they represent 
the simplest - or perhaps most pristine - example of a more general class of interplanetary 
coronal mass ejections (ICMEs). However, we emphasize that the theoretical relationship 
between ICMEs and MCs is not well understood. 

A simple, yet successful technique for modelin g flux ropes in the solar wind is the so- 
called force-free model (e.g., (Burlaga 1988; Lepping et al. 1990; Marubashi 1997; Mulligan 
et al. 2001; Hidalgo et al. 2002a)). In its early inception, clouds were envisaged to have 
cyUndricaUy-symmetric force-free configurations (e.g., Burlaga (1988)) and solutions of V x 
B = »B involving Bessel functions of the first kind (Lundquist 1950) yielded profiles for the 
magnetic field components that often matched observations well. Recent extensions to this 
have considered: (1) cloud expansion via ad hoc modifications to the model field components 
(e.g., Marubashi (1997); (2) multi-spacecraft observations (e.g., Mulligan et al. (2001); (3) 
including two radii of curvature - one for the flux-rope cross-section, and one for the curvature 
of the rope axis; (4) non-force-free effects by fitting to current densities (Hidalgo et al. 2002a) ; 
and (5) elliptical cross-sections (Hidalgo et al. 2002b; Russell & Mulligan 2002). 

An alternative approach for determining the properties of magnetic flux ropes has been 
developed recently by Hu & Sonnerup (2002). Using magnetic field and plasma measure- 
ments, they are able to estimate a number of parameters characterizing the flux rope, includ- 
ing it’s orientation, impact parameter, size, maximum field strength, and twist by integrating 
the non-linear, plane Grad-Shafranov equation. What is particularly appealing about this 
technique, as it pertains to the present study, is that it allows one to reconstruct the global 
structure of the cloud in the plane perpendicular to the cloud’s axis. So far, only 3 events 
have been analyzed, including a double flux rope (Hu & Sonnerup 2002; Hu et al. 2003). 

The success of the force-free models has led to a view of MCs as locally cylindrically- 
symmetric flux ropes that may or may not be connected back to the Sun. This is summa- 
rized by the illustrations shown in Figure 1 (a through c). Extensions to this picture have 
included the effects of solar rotation (Figure Id), and dynamical effects resulting from fast 
CMEs ploughing into ambient solar wind ahead and “flattening” (Figure le). Ulysses, while 
immersed in continuous fast, quiescent coronal-hole flow observed a distinct set of CME 
related signatures leading Gosling et al. (1994) to propose the picture shown in Figure 1(f). 

Fluid and MHD models have also been developed to explore the initiation (e.g., Mikic 
& Linker (1994); Linker &; Mikic (1997); Antiochos et al. (1999); Lin et al. (1998); Titov 
& Demoulin (1999); Filippov et al. (2001)) and evolution (e.g., Riley et al. (1997); Riley & 
Gosling (1998); Riley et al. (2002); Odstrcil et al. (2002); Vandas et al. (2002); Cargill & 
Schmidt (2002); Riley et al. (2003)) of CMEs near the Sun and in the solar wind. These 
models include a rich variety of physics and have been quite successful in reproducing a 
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wide range of observational signatures. However, as the level of sophistication in the model 
increases, so to does the difficulty in interpreting the results. In particular, it can be diffi cult 
to distinguish between competing processes. 

Previous studies have considered some aspects of kinematic distortion. Newkirk et al. 
(1981), for example, considered the temporal evolution of “bubble” -shaped structures, illus- 
trating primarily the azimuthal elongation of the bubble and distortion in response to fast 
and slow solar wind streams. Suess (1988) described the basic kinematic deformation of an 
initially circular cross section as it convected away from the Sun in a spherically-expanding 
solar wind. However, by adopting the prevalent view that the cross section of magnetic 
clouds was approximately circular, he reasoned that magnetic tension must be sufficiently 
strong to resist such a distortion. He further suggested that the observed radial expansion 
often seen in magnetic clouds could then be produced by the clouds tendency to minimize 
departures from force-free equilibrium brought about by spherical expansion. 

Multi-spacecraft observations of the same magnetic cloud offer an opportunity to infer a 
more global picture. Crooker &; Intriligator (1996), for example, analyzed a magnetic cloud 
observed by two spacecraft: IMP 8 at Earth and Pioneer 11 at 4.8 AU. The longitudinal 
separation of the two spacecraft was 30° and clearly indicated a distention of the ejecta. 
On the other hand, the observed declining speed profiles at the two spacecraft suggested 
that radial expansion had taken place at least close to the Sim. Thus they concluded that, 
in addition to kinematic distortion, dynamic expansion must have also taken place. Riley 
et al. (2003) studied a magnetic cloud observed by ACE at 1 AU and in the ecliptic plane 
and by Ulysses at 5 AU and S22° heliographic latitude. In conjunction with an MHD 
simulation, they derived a global picture of the event as it must have looked crossing the two 
spacecraft, consisting of a flattened, convex-outward ejecta where the transverse dimensions 
were significantly larger than the radial dimension. 

In this paper, we analyze the propagation and evolution of flux-rope CMEs using a 
kinematic approach. Our purpose is to show that, using the simplest technique possible, 
the most basic shape of a flux-rope CME in the solar wind is not cylindrical, but rather 
a convex-outward, “pancake” shape. Such a technique, while highly idealized, bridges the 
gap between the simple force-free models and the more complex numerical MHD models. 
Although we describe the evolution of the flux ropes in terms of their magnetic field, the 
results apply equally well to non-flux-rope CMEs. Our goals axe to isolate non-dynamical 
effects from the MHD solutions and describe the fundamental cross-sectional structure that 
we believe should be incorporated into future force-free, and non-force-free flux rope fitting 
models. The prevalent view that force-free flux ropes (and ICMEs in general) are cylindrical 
structures is unlikely to be met in most, if not all cases, implying that the illustrations 
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collected in Figure 1 are incorrect, even at the schematic level. 


2. Kinematic Evolution 

To model the kinematic evolution of magnetic flux ropes in the solar wind (where we 
use the term “solar wind” to refer to everything above the solar surface) we consider an 
initially-circular flux rope configuration, centered IRs above the solar surface, with an ini- 
tial radius of IRs (Figure 2a). Because of the implied symmetry in the azimuthal direction, 
by extension, this is a torus in three dimensions, and thus “field lines” ’ close back on them- 
selves and not back to the Sun. (Note that to be considered “cylindrical” would require that 
the radius of the flux rope, r/ iuirope <C IRs)- Although highly-idealized, numerous observa- 
tions of CMEs launched off the solar limb suggest that such a shape is reasonable starting 
point. We consider two types of kinematic evolution: (1) spherical expansion due to sim ple 
convection with the ambient solar wind in a diverging geometry; and (2) expansion due to a 
uniform pressure gradient between the flux rope and the ambient solar wind. The first effect 
is a true kinematic process, whereas the second is a dynamic process that we are treating 
kinematically. These two effects are shown in Figure 2a and b, respectively. Modeling the 
convective evolution involves nothing more than tracking the loci of a set of points in a uni- 
form spherical expansion. We assume purely radial flow so that the points maintain constant 
latitude, but move outward by a distance A R = v r x At, where A t is some arbitrary time step 
and v r is the velocity of the solar wind. Expansion due to a pressure gradient between the 
CME and the ambient solar wind can be implemented kinematically by calculating the unit 
outward normal to the CME boundary (i.e., the direction of —VP, ass uming the pressure of 
the flux rope is higher than that of the surrounding medium so that it expands rather than 
collapses) and shifting it by some constant value (i.e., AR = v exp At — (^)t^ xp x At x ex), 
where is some arbitrary constant expansion speed, n is the number of concentric circles 
that make up the flux rope, and i is the index of the circle under consideration. Thus the 
factor (jj) reflects the fact that pressure-gradient expansion is smallest near the center of 
the flux rope and largest at the edge. Note that a circle will maintain its shape under this 
operation, however, a more complex structure will deform. 

Practically, these ideas are easily incorporated into a simple computer program using 
a combination of translations, numerical derivatives, and transformations between cartesian 
and polar coordinates. We apply these operations to a set of concentric circles (as illustrated 
in Figure 2) to mimic the evolution of the poloidal component of the magnetic field wit hin 
the flux rope. 
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3. Results 

We consider first the effects of spherical expansion (i.e., convective evolution) alone on 
the flux rope. Figure 3 summarizes the evolution of an initially cylindrical flux rope at 
various heliocentric distances. Initially, the effect of the diverging spherical geometry is to 
flatten the trailing edge of the flux rope. Later on the flux rope develops a convex-outward, 
“pancake” shape. This is simply the result of the spherical geometry of the system and is not 
related to any interaction between the flux rope and the ambient wind (since there is none) . 
The initial radial width of the flux rope was IRs, and it remains this value throughout its 
evolution. By 1 AU (215 Rs) the flux rope’s transverse dimension dominates its radial width 
and it has the appearance of being “smeared” over a circular arc. The angular span of the 
flux rope remains a constant 60°. 

Observed MCs have a typical radial width of 0.23AU at 1 AU (Burlaga 1988) and so 
spherical expansion alone cannot produce such dimensions given their initial radial extent 
at the Sun. At the Sun, white light observations suggest that CMEs are undergoing strong 
expansion due to a higher internal pressure. In the solar wind, in situ observations often show 
declining speed profiles within ICMEs demonstrating that they are continuing to expand. 
However, where this expansion is most significant is not well known. Thus for simplicity, in 
our analysis we assume a linear expansion as a function of time. As we will see, comparison 
with MHD results suggests that this is a reasonable approximation to make. It is important 
to note that while successive operations of spherical expansion axe commutative, mixing 
pressure- gradient expansion with spherical expansion is not. This is a consequence of the 
fact that the direction of the pressure-gradient expansion is along the normal to the boundary 
of the flux rope, which is a function of time. Thus we evolve the flux rope by performing a 
repeated sequence of spherical expansion followed by pressure-gradient expansion. 

Figure 4 illustrates these effects on our idealized flux rope. We chose the pressure- 
gradient expansion velocity to be such that a radial width of approximately ~ 0.25 AU was 
produced at 1 AU. As in the previous case, the initial response of the flux rope to the 
spherical expansion is seen as a flattening of the tr ailing edge. At 3 Rs the aspect ratio 
(which we define to be the ratio of the latitudinal extent of the flux rope to its radial extent 
at the equator) is 1.4. By 30 Rs, however, the effects of pressure-gradient expansion can be 
clearly seen in increased radial width of the flux rope (~ 8 Rs) and an aspect ratio of 5.1. 
At 108i?s the aspect ratio has risen to 6.7. By 1 AU, the flux rope extends ~ 3b Rs (~ 0.25 
AU) in the radial direction and ~ 200i?s out of the equatorial plane, leading to an aspect 
ratio of 5.1. We can also look at the angular span of the ejecta. For the 5 snapshots shown 
in Figure 4, from top to bottom, they were: 66°, 98°, 114°, and 126°, respectively. T his 
variation is due entirely to the pressure-gradient expansion. 
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4. Comparison with an MHD Simulation 

We have developed a global MHD model of CME eruption at the Sun and evolution 
in the solar wind out to 5 AU (Odstrcil et al. 2002). Although idealized, we have found 
that the simulation reproduces many generic features of magnetic clouds (Riley et al. 2003) 
as well as predicting the presence of new phenomena that can be identified in in situ data 
(Riley et al. 2002). In Figure 5 we summarize the evolution of the flux rope at various times 
following its eruption. The radial lines indicate the latitudinal extent of the flux rope. A 
detailed description of these results has been presented elsewhere (Riley et al. 2002, 2003) 
and we restrict ourselves here to a brief description of the latitudinal and radial evolution 
of the flux rope. Comparison of these profiles with the kinematic evolution summarized in 
Figure 3 suggests that the two effects of spherical expansion and pressure-gradient expansion 
dominate the evolution of the large-scale structure of the flux rope. The major difference 
between the two profiles is due to the presence of the heliospheric plasma sheet, a region of 
slower, denser plasma in this case offset from the equatorial plane by —10°, which has the 
effect of squeezing the ejecta. This highlights that in reality, the presence of a more highly 
structured ambient solar wind will play a significant role in distorting the flux rope beyond 
the concepts discussed here. Defining the aspect ratio to be the ratio of the maximum 
vertical extent of the ejecta to its radial extent in the equatorial plane, we find it to be: 1.0, 
3.1, 5.8, and 9.3 for the 4 snapshots shown in Figure 5. Obviously these values are quite 
sensitive to the particular latitude chosen. Nevertheless, the variation from one time to the 
next is not. Thus the trend toward higher aspect ratios with increasing heliocentric distance 
is a real effect. We can also compute the angular span of the ejecta. Since the flux rope was 
not ejected exactly along the equatorial plane, but displaced southward by ~ 10° or so, we 
have computed the angle: 9 = 2 x tan -1 (A?// Ax), where Ay is the vertical span of the flux 
rope, and Ax is the horizontal span of the ejecta. These values are: 28°, 47°, 56°, and 59° 
for the 4 snapshots. Again, these numbers are not meant to be used for detailed quantitative 
analysis: It is their variation that primarily concerns us. We note that there is substantial 
angular expansion between 20 hours and 51 hours, when the flux rope moved from ~ 4.7 Rs 
to ~ 0.35 AU. This expansion had slowed by 99 horns (~ 0.9 AU), and still further by 435 
hours (~ A.AAU). 

We can apply force-free model fits to these simulation results by generating time series 
profiles at specific locations to mimic in situ measurements. In Figure 6 we compare such 
profiles with a force-free fit to the model data. The top 3 panels show speed, number density 
(n p ), and temperature (T p ) as a function of time. The hypothetical spacecraft was located 
at 2 °N heliographic latitude and 5 AU from the Sun. The boundaries of the flux rope are 
indicated by the vertical lines. The declining speed profile within the flux rope indicates that 
the structure is expanding. This is supported by the low density and temperature within the 


-7- 


ejecta. It is also driving a fast forward shock, st anding almost 1 day ahead of the leading 
edge of the flux rope. The remaining 4 panels summarize the magnetic field components as 
well as the magnitude of the field. They are shown in the rtn coordinate system where e r is 
the radial direction (positive is away from the Sun), e< is parallel to the equatorial plane and 
in the direction of planetary motion, and e„ completes the right-handed system. (Relative to 
the more familiar spherical coordinate system ( r,6,<j> ), e t = and e„ = —eg). The simulated 

magnetic components show the classic features of a magnetic cloud with its axis lying parallel 
to the equatorial plane and perpendicular to the radial direction: The radial component of 
the field remains zero throughout; the azimuthal component rotates from zero through to a 
maximum and back to zero, while the meridional component falls from a maximum , through 
zero, and to a minim u m negative value. The azimuthal component reflects changes in the 
axial component of the field and the meridional component reflects changes in the poloidal 
component. The force-free fitting technique we have applied follows that of Lepping et al. 
(1990) with the additional simplification that we know the precise orientation of the flux 
rope and so do not have to include these parameters in the fit. The magnitudes of the field 
components are set by the maximum value of the observed (or simulated in this case) a xial 
field. Qualitatively, there is a fair agreement between the s imula ted profiles and the Bessel 
function fits. Of notable exception are: (1) the asymmetry present in the simulated profiles; 
(2) the larger extrema in the simulated B n profiles; and (3) the relatively flat (but slightly 
declining) simulated magnetic field magnitude profile. The asymmetry is due primarily to 
the expansion of the ejecta, while the larger simulated B n profile and flat/declining field 
magnitude profile are due to kinematic distortion. 


5. Summary and Discussion 

In this study, we have presented a kinematic study of the evolution of CMEs in the solar 
wind. By isolating the effects of spherical expansion and pressure-gradient expansion, and 
comparing the results with MHD simulations, we have been able to illustrate the importance 
of kinematic evolution on CME morphology. 

One of the most important dynamical effects included in MHD simulations, and not 
present in our kinematic analysis, is due to relative speed difference between the flux rope 
and the surrounding medium. CMEs are often launched with speeds significantly faster 
than the ambient solar wind, driving a forward wave, which eventually steepens into a shock 
and generates a region of enhanced pressure behind it (the sheath). The pressure gradients 
associated with this act to flatten the flux rope profile further on its leading edge. In contrast, 
in the region behind the flux rope an expansion wave (rarefaction region) forms, accelerating 
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the slower wind behind the flux rope and decelerating its trailing edge. This acts to resist 
the development of the convex-outward shape at the trailing edge. A second important effect 

- and one that we cannot address within the confines of the kinetic approach described here 

- concerns the properties of the ambient solar wind. At solar minimum , for example, when a 
band of slow, dense wind emanates from lower latitudes, while fast tenuous wind flows from 
large polar coronal holes, a CME, initially launched into the slower wind will likely penetrate 
into the faster flow. The evolution of the flux rope in these two different environments will 
be fundamentally different and the resultant global morphology of the flux rope will be 
significantly more complicated than suggested here (e.g., Riley et al. (1997)). Magnetic 
reconnection at regions of strong velocity shear could further complicate the evolution of the 
CME by altering its initial topology (Schmidt & Cargill 2001). 

Our kinematic examples show that when only spherical expansion is considered, the flux 
rope maintains constant latitude (~ 30° in the case shown in Figure 2). On the other hand, 
when pressure-gradient expansion is added, the flux rope penetrates beyond ~ 60° latitude. 
In reality, the expansion associated with pressure gradients is unlikely to be a simple linear 
function of time as we have assumed. In fact, from Figure 5 we infer that from panels 1 
through 4 the angular expansion progressed from 70% to 18% to 6%. Since the latitu dinal 
extent of the flux rope is very sensitive to where this expansion takes place (due to the 
diverging geometry of the system), the true latitudinal extent may be substantially different, 
and our kinematic results should be taken as indicative, but not quantitatively correct. 

In our kinematic analysis we have ignored any role played by magnetic tension in re- 
sisting the distortion of the cloud by spherical expansion (Suess 1988). This force, however, 
is included in the MHD model. Thus the fact that the MHD simulations show the same 
qualitative distortions as the kinematic treatment suggests that magnetic tension does not 
provide a dominant restoring force. Moreover, for ICMEs that (1) do not contain well-defined 
flux ropes; and/or (2) are not low-beta structures, it is even more unlikely that tension plays 
a significant role in resisting distortion due to spherical expansion. It has also been shown 
that magnetic tension provides very little resistance against strong velocity shear (Schmidt 
& Cargill 2001). Force free models can in principle deduce the contribution played by mag- 
netic tension. The differences between the fitted force-free profiles and the observed profiles 
provides a measure of the non-circularity of the poloidal component of the magnetic field. 
The asymmetry in the field magnitude in Figure 6 is, in a large part, due to the flattening 
of the ejecta. The effect is also seen in the meridional field component (B n ), where the 
force-free model underestimates the field. The flattening effect of spherical expansion adds 
to the meridional component, without affecting the maximum value of the axial field. 

These kinematic distortions are a basic feature of all MCs. Cylindrical force-free models, 
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including those which treat the effects of expansion along the trajectory of the spacecraft, 
may benefit greatly by incorporating the effects of the evolution described here. One way to 
accomplish this would be to add an “aspect ratio” free parameter into the fitting technique. 
This might be defined as the ratio of the semi-major axis to semi- min or axis of the flux rope 
(in a coordinate system based on an arc sweeping in latitude through the flux rope). Since the 
effect of spherical expansion is to increase the meridional (or 6 ) component of the magnetic 
field at the expense of the radial component, this new free parameter could be incorporated 
here. Models of flux ropes that fit to an elliptical geometry (e.g., Hidalgo et al. (2002b)), 
while an improvement over cylindrical fits, do not take into account the curvature of the 
MC in the meridional plane. In fact, to generate a truly elliptical structure at 1 AU would 
require that the northern and southern flanks of the ejecta were traveling faster than lower 
latitude portions, at just the right amount to maintain the elliptical geometry - a scenario 
that is unlikely to be met in any typical ambient solar wind. In spite of these limitations, 
force-free and non-force-free models: (1) are easy to apply to in situ measurements; and (2) 
provide the only way to infer the cloud’s orientation, chirality, and radial dim ension. Thus 
they will likely remain an important tool for analyzing magnetic clouds for the foreseeable 
future. It should be emphasized, however, such techniques apply only locally in the vicinity 
of the cloud. In the presence of a structured ambient solar wind, the global picture of the 
ejecta will likely be significantly different from that drawn from the fitted parameters. The 
accuracy of these fitting techniques, both in view of the results presented here as well as 
from dynamic contributions, remains to be established. We have recently begun a study 
using simulated time series (extracted from global MHD simulations with realistic ambient 
solar wind conditions), where we know the local and global properties of the flux rope, to 
quantify the accuracy of a variety of these fitting techniques. 

In closing, we reiterate that while the techniques employed in this study axe exceptionally 
simple, yet their implications for MC and ICME modeling are important. On one hand, 
they compel future force-free methods to include the convex-outward geometries that are 
an inescapable consequence of the spherical expansion of the CME and the solar wind in 
which it is embedded. And on the other hand, they hi ghlight the pronounced contribution 
of kinematic effects to CME evolution that have often been ascribed to dynamical effects. 
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Fig. 1. — An illustrative collection of flux rope schematics, highlighting different aspects of 
the properties and evolution of ICMEs. (a) Global profile of a magnetic cloud approaching 1 
AU (Burlaga et al. 1990); (b) Helical properties of cylindrical flux rope field lines (Bothmer 
& Schwenn 1998); (c) Spacecraft trajectory through a cylindrical flux rope (Lepping et al. 
1990); (d) Global picture of flux rope approaching 1 AU, including the effects of solar rotation 
(Marubashi 1997); (e) Dynamical effects on evolution of a fast CME (Gosling 1990); and (f) 
Schematic of CME evolution at high heliographic latitudes (Gosling et al. 1994). 

Fig. 2. — Kinematic effects on evolution of magnetic flux rope, (a) Convective evolution 
or spherical expansion; and (b) expansion of the flux rope caused by a pressure gradient 
between it and the ambient solar wind. 

Fig. 3. — Evolution of a flux rope due only to the effects of spherical expansion. Shown are 
snapshots of the flux rope at ~ 3, 32, 110, and 215 Rs- 

Fig. 4. — Evolution of a flux rope due to the effects of both spherical expansion and pressure- 
gradient expansion. Shown are snapshots of the flux rope at ~ 3, 32, 110, and 215 Rs- 

Fig. 5. — Radial velocity (color contours), magnetic flux function (black, closed contours), 
and number density (red contours) for 3 times during an MHD calculation. 

Fig. 6. — Comparison of MHD results with force-free fitting to magnetic field parameters for 
the CME shown in Figure 5. 
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